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Abstract 

We present a simple method for determining the (correlated) uncertainties of 
the light element abundances expected from big bang nucleosynthesis, which 
avoids the need for lengthy Monte Carlo simulations. Our approach helps to 
clarify the role of the different nuclear reactions contributing to a particu- 
lar elemental abundance and makes it easy to implement energy-independent 
changes in the measured reaction rates. As an application, we demonstrate 
how this method simplifies the statistical estimation of the nucleon-to-photon 
ratio through comparison of the standard BBN predictions with the observa- 
tionally inferred abundances. 

26.35. +c, 98.80.Ft 
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I. INTRODUCTION 



Big bang nucleosynthesis is entering the precision era |IJ. On the one hand, there has been 
major progress in the observational determination of the abundances of the light elements 
D IIH, 3 He jU|, 4 He |||7|], and 7 Li ||||, although the increasing precision has highlighted 
discrepancies between different measurements (see Refs. HlOHT^] for recent assessments). 
Secondly, we have a sound analytical understanding of the physical processes involved |T3|JT^| 
and the standard BBN computer code [|15|,|nj which incorporates this physics is robust and 
can be easily altered to accomodate changes in the input parameters, e.g. nuclear reaction 
rates [JD]]. The comparison of increasingly accurate observationally inferred and theoretical 
abundances will further constrain the values of fundamental parameters, such as the nucleon 
density parameter (see, e.g., Ref. JTSI) or extra degrees of freedom related to possible new 



physics beyond the Standard Model (see, e.g., Ref. |1|). It goes without saying that error 
evaluation represents an essential part of such comparisons. 

Because of the complex interplay between different nuclear reactions, it is not straight- 
forward to assess the effect on a particular elemental yield of the uncertainties in the exper- 
imentally determined reaction rates. The authors of Ref. [EJ first employed Monte Carlo 
methods to sample the error distributions of the relevant reaction cross-sections which were 
then used as inputs to the standard BBN computer code. This enables well-defined confi- 
dence levels to be attached to the theoretically predicted abundances, e.g. the abundance 
range within which say 95% of the computed values fall correspond to 95% CL. limits on the 
expected abundance. It was later realized that error correlations are also relevant, and can 
be estimated with the same technique [pTjp2| . The Monte Carlo (MC) approach has since 
become the standard tool for comparing theory and data |iTp3H 2l| . However, although it 
can include refinements such as asymmetric or temperature-dependent reaction rate uncer- 
tainties |T7[ , it requires lengthy calculations which need to be repeated each time (any of) 
the input parameters are changed or updated. Since we may expect continued improvement 
in the determination of the relevant parameters, it is desirable to have a faster method for 
error evaluation and comparison with observations. 

In this work we propose a simple method for estimation of the BBN abundance un- 
certainties and their correlations which requires little computational effort. The method, 
based on linear error propagation, is described in Sec. II. A concrete application is given in 
Sec. Ill, where theory and observations are compared using simple \ 2 statistics to obtain 
the best-fit value of the nucleon-to-photon ratio. In Sec. IV we study with this method the 
relative importance of different nuclear reactions in determining the synthesized abundances. 
Conclusions and perspectives for further work are presented in Sec. V. 



II. PROPAGATING INPUT CROSS SECTION UNCERTAINTIES TO OUTPUT 

ELEMENTAL ABUNDANCES 

A. Notation and input 

The four relevant element abundances Yj considered in this work are defined in Table I. 
(Note that the abundance of 4 He is conventionally quoted as a mass fraction, while the 
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abundances of D, 3 He and 7 Li are ratios by number.) In BBN calculations, the Yj's depend 
both on model parameters (the nucleon-to-photon ratio T], the number of neutrino families 
and on a network of nuclear reactions Rj.: 



N v , etc. 



Yi = Yi(ri, N v 



{Rk}) 



(i) 



The most important R^s are listed, numbered as in Ref. jig] , in the first two columns 
of Table II, while our default inputs for the rates Rk and their la uncertainties ±A Rk 
are given in the third and fourth columns. The numerical values are given in ratio to the 
reference reaction rates compiled in Table 1 of Ref. |17]; we have chosen identical values 
(i.e., Rk = 1) except for Ri, the neutron decay rate, where we adopt the most recent world 
average for the neutron lifetime of r n 
888.54 ± 3.73 s used in Ref. 



886.7 ± 1.9 s 27,28 



T, 



T7T 



as compared to the value of 
The fractional uncertainties ±ARk/Rk (k ^ 1) have 
also been taken from Ref. |l7j (see their Table 2), assuming conservatively the largest value 
for the temperature-dependent errors AR 7 and Ai? 10 .Q 



B. The method 

For simplicity we consider only the standard BBN case (i.e., N v = 3, etc.), so that r\ is 
the only model parameter being varied in the calculation of the abundances Yi and of their 
uncertainties af. 

Y i = Y i (Tj)±<T i (T l ) - (2) 

Our method can however be easily generalized to nonstandard cases. 

For a relatively small change 5 Rk of the input rate Rk {Rk - > Rk + SR k ), the corre- 
sponding deviation SYi of the z-th elemental abundance (Yi — > Yi + 5Yi), as given by linear 
propagation, reads 

SYi( V )=Yi( V )^:X ik ( V ) 5 -^, (3) 

k R k 

where the functions \ik(ff) represent the logarithmic derivatives of Yi with respect to Rk'- 

, , . dlnYAri) 

In general, the deviations SYi in Eq. (^) are correlated, since they all originate from 
the same set of reaction rate shifts {S Rk}- The global information is contained in the error 
matrix (also called covariance matrix) [p9fl ) which is a generalization of the "error vector" S Yi 
in Eq. (|3]). In particular, the abundance error matrix crfArj) obtained by linearly propagating 
the input ±ler reaction rate uncertainties ±AR k to the output abundances Y$ reads: 



1 0ur method for error propagation requires that the AR^/R^s be constant (i.e., temperature 
independent). We will comment on this point at the end of this Section. 
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= Yi{v)YMEMv)^k(v) (^) 2 • (5) 

This matrix completely defines the abundance uncertainties. In particular, the la abundance 
errors cjj of Eq. (0) are given by the square roots of the diagonal elements, 



°i(v) = \/ vfM , (6) 
while the error correlations pij can be derived from Eqs. (|)]||) through the standard definition 

„..(„)= ^ (7) 

Thus Eqs. represent all that is required to calculate the errors in the predicted abun- 
dances and their correlations. 

Note that the relevant physics is contained entirely in the central values Y,. and in their 
logarithmic derivatives A^, which have to be evaluated just once with a BBN numerical 
code, thus dramatically reducing the required computing time.f] We have made a further 
check of the linearity of the error propagation by calculating the logarithmic derivatives with 
increments equal to ARk (default) and 2 ARk, obtaining practically the same functions \ik 
in either case. This means that doubling the error on Rk also doubles the corresponding 
error component of Y i: i.e. the error propagation is indeed linear. 

We think it useful to present the results of this exercise in the form of tables so all 
calculations that will now follow can be done on a pocket calculator. Tables III and IV show 
the coefficients of polynomial fits to Y{ and \k, respectively, for rj in the usually considered 
range 10~ 10 — 10~ 9 .[] The abundances Y^rf) with their associated ±2 standard deviation 
error bands calculated through Eq. (|6]) are shown in Fig. [I]. The functions \ik(v) are shown 
in Fig. ^; note that some of these vary strongly (and even change sign) with r], indicating 
that the physical dependence of the Y^'s on the Rk's may be quite subtle. Although some 
general features of this dependence have been addressed in Ref . [p| , further work is needed to 
interpret the functional form of the Ajfc's in Fig. |2|. We intend to address this issue elsewhere. 
Finally, Fig. ^displays the fractional uncertainties Oi/Y* and their correlations as derived 
from Eqs. (|5^|7|). Notice that, in general, the error correlations are non-negligible and should 



be properly taken into account in statistical analyses, as first emphasized in Ref. [21 . 



2 The logarithmic derivatives are numerically defined as A^(r/) = \Yi(r), Rk + ARk) — Yi(r], Rk — 
A Rk)]Rk/2Yi(r], Rk)A Rk, at given r]. We find negligible difference between left and right 
derivatives. 

3 Small corrections to the helium abundance Y4 due to Coulomb, radiative and finite temperature 
effects, finite nucleon mass effects and differential neutrino heating, have been incorporated accord- 
ing to the prescription given in Ref. (T^[ . We have used the BBN code [|jl|] with the lowest possible 
settings of the time steps in the (2nd order) Runge-Kutta routine, which allows rapid convergence 
to within 0.01% of the true value JkT]. We understand that our results are in good agreement with 



a recent independent computation of Y4 using a new BBN computer code [31|. 
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In summary, the recipe for evaluating the BBN uncertainties ±<7j affecting the y's for a 
given value of rj is: 

(z) Determine the abundances Yiirf) and their logarithmic derivatives Xik(v) using Tables 
III and IV, respectively; 

(ii) If the central values of the reaction rates Rk are updated (Rk — > Rk+8 Rk) with respect 
to those reported in Table II, then update also the central values of the abundances 
(y.^y. + 5 y.) through Eq. ©; 

(in) For given reaction rate uncertainties AR^. (e.g., from Table III), compute the abun- 
dance errors cr, and their correlations p^- using Eqs. (0-0). 

C. Comparison with MC estimates and remarks 

Our approach is based on the linear propagation of errors originating from many inde- 
pendent sources (i.e. the Rk's). One can expect that this method will work reasonably well, 
both because the input fractional uncertainties AR^/Rk are relatively small, and because 
the final output uncertainties <7j affecting the abundances y are "regularized" by the central 
limit theorem. Indeed, our ±2cr bands in Fig. [I] compare well with the MC-estimated bands 



of Refs. |17| , |21] , |23| -Pq1 , with small relative differences which depend, in part on different in- 
put Rk ± Ai?fc's, and that are not larger than the spread among the various MC estimates 
themselves. 

In order to be more quantitative, we compare in Fig. ^ (upper panel) the MC evaluation 
of the fractional uncertainties <Ji/Yi as derived from Ref. |T7| P| with our analytic estimate 
(using, for this exercise, the same input parameters). There is good agreement between 
these two totally independent estimates. In Fig. |4] (lower panel) we also show a comparison 
with the only MC evaluation of we are aware of (viz., Ref. p2|), obtaining again good 
agreement with our calculation when the same input Rk ± ARk are used. We conclude the 
discussion of Fig. [| by noting that the uncertainty and the correlations P(2+3)j related 



to the often-used combination of abundances Yi 2 +3) = V 2 + Y 3 = (D + 3 He)/H are given, 



within our approach, by: 



a 



J2 _2 i 2 



(2+3) 



a 2 + a 3 + 2p 23 a 2 a 3 , (8) 

P(2+3)i 0r (2+3)° r j = P2j Cr 2 CT j +P3j a 3 a j ■ ( 9 ) 

There are, of course, some refined features of the MC approach that cannot be addressed 
with our method, such as asymmetric or temperature-dependent uncertainties A Rk [17 . 
However we consider these refinements not essential for practical applications. In a sense, 
the possible asymmetry between "upper" and "lower" errors is where one wants it to be. For 



4 The MC values of cji/Yi have been read off the (small) panels of Fig. 27 in Ref. [17| and, therefore, 
may be subject to small transcription errors. 
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instance, if one assumes a priori symmetric errors in the astrophysical S'-factors, then asym- 
metric errors are induced in the thermally averaged reaction rates R ~ (av) ; conversely, the 
requirement of a priori symmetric A errors requires that the input 5-factor uncertainties 
are readjusted, as discussed in Ref. Ul7j . Although the authors of Ref. |17]] have adopted the 



latter option (| + ARk\ = | — Air^]), the former option or others are equally acceptable, 
and would clearly produce different outputs for the MC estimate of the abundance errors. 
For instance, the upper and lower errors of Y? appear to be rather symmetrical in the MC 



calculation of Ref. [17], while they are noticeably asymmetrical in Ref. 



Concerning the temperature-dependent |nj uncertainties AR7 and Ai? 10 , which affect 
mainly the estimate of Y 7 , our conservative choice in Table II proves to be successful for the 
estimate of o-jjY-j (see Fig. £|). This seems to indicate that the uncertainties at low temper- 
atures (which are larger than those at high temperatures (T^j) dominate in the estimate of 
these errors. In any case, since practically all reaction rate uncertainties A R^ contribute to 
the final value of a 7 , temperature-dependent refinements in the propagation of just two of 
these (A R 7 and A R 10 ) do not appear to be decisive for the estimate of the global error a 7 . 

In conclusion, we have shown that our simple analytic method for error evaluation repre- 
sents an useful alternative to lengthy and computationally expensive MC simulations. Both 
the magnitude and the correlations of the total errors affecting the theoretical abundances 
are reproduced with good accuracy. We therefore advocate the use of this method for BBN 
analyses as an alternative to MC simulations.^ 



III. DETERMINING THE LIKELY NUCLEON-TO-PHOTON RATIO 

The comparison of the predicted primordial abundances Y^rj) ±<jj(?7) with their observa- 
tionally inferred values Y \ ± ex, through a statistical test allows extraction of the likelihood 
range for the fundamental parameter r]. So far, this has been done either through fit-by-eye 
(see, e.g., Ref. fl7| ) or by Monte Carlo-based maximum likelihood methods (see, e.g., Refs. 
26]). In this Section we show how limits on r] can be simply extracted using \ 2 statistics 



based on the method described in the previous Section. 

Assuming that the errors Zfj in the determinations of different abundances Y i are uncor- 
rected, the experimental squared error matrix of- is simply 

0% = SijWiWj , (10) 

where 8^ is Kronecker's delta. The total (experimental + theoretical) error matrix Sfj is 
then obtained by summing the matrices in Eqs. ([5|JT0|): 

Syfo) = 4(^+4- (11) 



5 A parallel situation holds in the field of solar neutrino physics, where the correlated uncertainties 
of the neutrino fluxes predicted by solar models have been estimated through both Monte Carlo 
simulations [32] and linear propagation of input errors [33]. The latter technique has proved more 
popular because of its ease of use. 
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Its inverse defines the weight matrix Wij(rf): 

= [sUvT 1 - (12) 

The x 2 statistic associated with the difference between theoretical (Yj) and observational 
(Yi) light element abundance determinations is then [pSfl : 



X 2 (V) = EK(^) " Y i\ ■ W M ■ \Vj(v) - Yj\ . (13) 

Minimization of the x 2 gives the most probable value of i], while the intervals defined by 
X 2 — Xmin + A x 2 gi ye t ne likely ranges of rj at the confidence level set by A x 2 (for one 
degree of freedom, 77) .0 

In order to illustrate this, we estimate rj using recent observational data for the three light 
element abundances (Y 2 , 14, 17) whose primordial origin is most secure. It is well known that 
the observationally inferred values Y 2 and Y4 are still controversial, and the conflict between 
different determinations has driven a lively debate on the status of BBN (see Refs. |18|J34]| 
and references therein). In this paper we do not enter into this debate but rather apply our 
method to two possible (although mutually incompatible) selections of measurements which 
we name data set "A" and data set "B" : 

= 1.9 ±0.4 x 10~ 4 , 
Data set A : {% = 0.234 ± 0.0054 , (14) 
= 1.6 ±0.36 x 10~ 10 ; 



Data set B : 




3.40 ±0.25 x 10~ 5 , 
0.243 ± 0.003 , 
1.73 ±0.12 x 10~ 10 



(15) 



The data set "A" is used in Ref. [pqj , the authors of which make a detailed MC-based 
fit to rj, thus enabling comparison with our method. Note that their adopted value of 
the primordial deuterium abundance from observations of high redshift quasar absorption 
systems is consistent with another recent observation, Y 2 = (2.15 ± 0.35) x 10~ 4 M, but in 
conflict with the significantly smaller value reported in Ref. 0] , which we adopt for data set 
"B" . Similarly the primordial helium mass fraction inferred from observations of metal-poor 
blue compact galaxies which we adopt for data set "B" is from Ref. in which it is argued 
that previous analyses leading to the smaller value of Y 2 used in data set 'A" underestimate 
the true abundance (although this is disputed in Ref. 0.) Finally the estimates for the 
primordial lithium abundance in both data sets are based on observations of pop II stars, 
with the slightly higher value of Y 7 in data set "B" taken from an updated analysis. 
Readers who prefer to adopt different combinations of these, or indeed other, estimates for 



6 We remind the reader that A% 2 = 1, 2.71, 3.84, and 6.64 correspond to confidence levels (C.L.'s) 
of 68%, 90%, 95%, and 99%, respectively. 
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the primordial abundances are invited to perform their own fit to 77 by following the simple 
prescription given here. 

Before performing the x 2 fit, it is useful to get an idea of what one should expect by 
comparing the data with the theoretical predictions at various values of 77. Fig. ^| shows 
the theoretical predictions for the abundances F 2 , I4, and Y 7 in the three possible planes 
(Yi, Yj), for representative values of x = log 10 (77/10~ 10 ). The corresponding la error ellipses 
show clearly the size and the correlation of the "theoretical" errors. The observational data 
sets "A" and "B" are also indicated on the figure, as crosses with la error bars. Clearly the 
former prefers rj ~ 2 x 10~ 10 while the latter favours 77 ~ (4-5) x 10~ 10 . 

A more precise estimate of the likely range of 77 is obtained, as anticipated, through 
a x 2 fit- The results are shown in Fig. || The value of xLm * s almost zero for the fit 
to data set "A", indicating very good agreement between theory and observations, while 
it is somewhat larger for data set "B". Note that the characteristic minimum in the 7 Li 
curve at 77 rs 2.6 x 1CT 10 (see bottom panel of Fig. [I]) allows its measured abundance to be 
compatible with both high D/low 4 He (data set "A") or low D/high 4 He (data set "B"). 
The 95% CL. ranges allowed by each of the two data sets, as obtained by cutting the curves 
at Ax 2 = X 2 ~ Xmin = 3.84, are: 

Data set A : 77 = 1.78±gj 4 x KT 10 , (16) 
Data set B : 77 = 5.13±g;^ x KT 10 . (17) 

The range for case "A" agrees very well with the 95% CL. range estimated in Ref. f26|| 
with the same inputs but with a different method (Monte Carlo + maximum likelihood). Of 
course, the incompatibility between the above two ranges of 77 reflects the incompatibility 
between the input abundance data within their stated errors. 



IV. ROLE OF DIFFERENT REACTIONS IN LIGHT ELEMENT 

NUCLEOSYNTHESIS 

The role of the different nuclear reactions rates listed in Table II in the synthesis of 
the light elements can be studied by "perturbing" the values of the input reaction rates 
and observing their effect on the predicted abundances. More precisely, one can study the 
contribution to the total uncertainty <7, of Yj induced by a +la shift of R k ~- 

R k ^R k + AR k =► Yi -> Y + 5Y . (18) 

Within our approach, this can be done very easily using Eq. (3), with the AR k s from 
Table II. Of course, the results depend on the value chosen for 77. To illustrate various 
trends, we choose the best-fit values 77 = 1.78 x 10~ 10 and 77 = 5.13 x 10 -10 , corresponding 
to data sets "A" and "B" respectively. 

Figs. ^ and [5] show the deviations 5Yi (normalized to the total error o~j) induced by 
+lcr shifts in the R k s, plotted in the same set of planes as used for Fig. ||. The la error 
ellipses shown in these figures are obtained by combining the deviation vectors SYi/ai in 
an uncorrelated manner. Several interesting conclusions can be drawn from this exercise. 
As expected, the uncertainty in the weak interaction rate Ri has the greatest impact on 
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Y4 for the high value of r\ (Fig. §), since essentially all neutrons end up being bound in 
4 He. However at the lower value of 77 (Fig. |[|), the uncertainty in R 2 — the "deuterium 
bottleneck" - - plays an equally important role as R\ in determining F4 because nuclear 
burning is less complete here than at high 77. Similarly with reference to the reaction rates 
i?7, R\o — R12 which synthesize 7 Li, at low rj it is the competition between R~j and Rn 
which largely determines Y?, while at high 77 it is the competition between R\q and R12. The 
anticorrelation between I4 and Y% is driven mainly by R2 at low rj and, to a lesser extent, 
by -R4 and -R5, while the reverse is the case at high rj. The anticorrelation between Y4 and 
Y-j at low 77 is also basically driven by R2, while the correlation at high 7/ is due to both 
i?2 and -R4. Thus we have a direct visual basis for assessing in what direction the output 
abundances Yi are pulled by possible changes in the input cross sections Rk- 



V. CONCLUSIONS 

We have shown that a simple method based on linear error propagation allows us to 
quantify the uncertainties associated with the elemental abundances expected from big bang 
nucleosynthesis, in excellent agreement with the results obtained from Monte Carlo simu- 
lations. This method makes transparent which nuclear reaction rate is mainly responsible 
for the uncertainty in the abundance of a given element. If determinations of the primor- 
dial abundances improve to the point where the observational errors become smaller than 
the theoretical uncertainties (say for 7 Li), this will enable attention to be focussed on the 
particular reaction rate whose value needs to be experimentally better known. 

We have also demonstrated that for standard BBN, our method enables the use of simple 
X 2 statistics to obtain the best-fit value of rj from the comparison of theory and observations. 
At present there are conflicting claims regarding the primordial abundances of, particularly, 
D and 4 He, and different choices of input data sets imply values of 7/ differing by a factor of 
~ 3. However this quantity can also be determined through measurements of the angular 
anisotropy of the cosmic microwave background (CMB) on small angular scales. Within a 
decade the forthcoming all-sky surveyors MAP and PLANCK are expected to pinpoint the 
nucleon density to within ~ 5% pgfl . Such measurements probe the acoustic oscillations 



of the coupled photon-matter plasma at the (re) combination epoch and will thus provide 
an independent check of BBN, assuming rj did not change significantly between the two 
epochs.^ Nevertheless precise measurements of light element abundances, particularly 4 He, 
are still crucial because they provide a unique probe of physical conditions, in particular 
the expansion rate at the BBN epoch. To illustrate, if rj was determined by the CMB 
measurements to be ~ 2 x 10~ 10 (consistent with data set 'A"), but the abundance of 4 He 
was established to be actually closer to its higher value of ~ 24% in data set "B", this 



7 New physics beyond the Standard Model can change 77, e.g. by increasing the photon number 
through massive particle decay [36| or, more exotically, by decreasing the photon number through 
photon mixing with a shadow sector |$7||. However such possibilities are strongly constrained by 
the absence of distortions in the Planck spectrum of the CMB |38) and also, in the latter case, by 
the absence of Sakharov oscillations in the power spectrum of large-scale structure [39|. 
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would be a strong indication that the expansion rate during BBN was higher than in the 
standard case with N v = 3 neutrinos. Although the number of SU{2) doublet neutrinos is 
indeed 3, there are many light particles expected in extensions of the Standard Model, e.g. 



singlet neutrinos, which can speed up the expansion rate during nucleosynthesis ||19|| . The 
generalization of our method to such non-standard cases is straightforward and we intend to 
present these results in a future publication [[|(|] . It is clear that BBN analyses will continue 
to be important in this regard for both particle physics and cosmology. 
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TABLES 



TABLE I. The four light elemental abundances Yi considered in this work. Alternative symbols 
used in the literature are indicated in parentheses. 



Symbol ... or Definition 

Y4 (Y p ) 4 He mass fraction 

Y 2 (y2 P ) D/H (by number) 

Y 3 (y 3p ) 3 He/H (by number) 

Y7 {yip) 7 Li/H (by number) 



TABLE II. The BBN reaction rates and their la uncertainties ±A adopted in this work. 
The numbering follows Ref. [16 while the reference "unit" values (Rk = 1) correspond to the rates 
in Ref. Q. 



k_ 
T~ 

2 

3 

4 

5 

6 

7 

8 

9 
10 
11 
12 



Reaction 



Ri 



±AR k 



n — ► p e u e 

p(n,7)d 
d(p, 7 ) 3 He 
d(d,n) 3 He 

d(d,p)t 
t(d,n) 4 He 
t(a,7) 7 Li 
3 He(n,p)t 
3 He(d,p) 4 He 
3 He(a, 7 ) 7 Be 
7 Li(p,a) 4 He 
7 Be(n,p) 7 Li 



0.9979 
1 
1 
1 
1 
1 
1 
1 
1 
1 
1 
1 



±0.0021 
±0.07 
±0.10 
±0.10 
±0.10 
±0.08 
±0.26 
±0.10 
±0.08 
±0.16 
±0.08 
±0.09 



TABLE III. Polynomial fit to the central value of the elemental abundances, 
Yi = ao + a\x + CL2X 2 + CI3X 3 + CI43; 4 + 05a; 5 , with x = log 10 (77/10 -10 ) in the range 0-1. The 
abundances were obtained using the BBN computer code |l6| with the input R^s as in Table II. 
The value of I4 has been corrected using the prescription of Ref. [19]. The accuracy of the fit is 
better than 1/25 of the total theoretical uncertainty for each Y{. 



Y 2 x 10 3 
y 3 x 10 5 
y 4 x 10 1 
Y 7 x 10 9 



+0.4808 
±3.4308 
±2.2305 
±0.5369 



ai 
-1.8112 
-6.1701 
±0.5479 
-2.8036 



02 

+3.2564 
+8.1311 
-0.6050 
+7.6983 



Q3 
-3.3525 
-9.7612 
+0.6261 
-12.571 



a 1 

+ 1.8834 
+7.7018 
-0.3713 
+ 12.085 



Q5 
-0.4458 
-2.5244 
+0.0949 
-3.8632 
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TABLE IV. Polynomial fits to the logarithmic derivatives, ^ik( x ) = 
ao + a±x + ci2X 2 + a^x 3 + a^x 4 + a^x 5 , as functions of x = log 10 (?7/10 _lcl ) € [0, 1] (see Fig. |2|). 
In most cases, polynomials of degree < 5 provide sufficently accurate fits. Only non-negligible 
logarithmic derivatives are tabulated. 





k 


ao 


ai 


a.2 


03 


(24 


Ct5 


^2k 


1 


+0.7130 


-0.7964 


+0.4577 


+0.0914 










2 


-0.7025 


+0.2611 


+1.2008 


-0.8934 










3 


-0.0189 


-0.1879 


+0.2502 


-0.6806 










4 


-0.4228 


-0.1698 


+0.0207 


+0.0247 










5 


-0.4138 


-0.1477 


+0.1103 


+0.0010 










8 


-0.0073 


-0.0003 


+0.0801 


-0.0416 










9 


-0.0011 


-0.0348 


+0.0592 


-0.0511 








^3k 


1 


+0.0940 


+0.1892 


-0.1484 


-0.0199 










2 


+0.0981 


+0.9948 


-3.1667 


+3.4108 


-1.2845 







3 


+0.0610 


+0.1640 


+0.5368 


-0.2605 










4 


+0.3050 


+0.0805 


-0.4208 


+0.1555 










5 


-0.5118 


+0.1274 


+0.4081 


-0.2106 










6 


-0.0327 


+0.0829 


-0.0939 


+0.0362 










8 


-0.5580 


-0.0287 


+1.3574 


-0.8735 










9 


-0.1080 


-0.5089 


-1.0157 


+0.8163 








Mk 


1 


+0.8138 


-0.1465 


+0.0408 













2 


+0.0610 


-0.1962 


+0.2416 


-0.1049 










4 


+0.0082 


-0.0058 


+0.0034 













5 


+0.0075 


-0.0058 


+0.0034 











^7k 


1 


+1.9638 


+1.8520 


-19.721 


+27.542 


-11.041 







2 


-0.9214 


-1.6472 


+5.6187 


+54.059 


-112.80 


+56.426 




3 


+0.0500 


-1.0433 


+4.1384 


-2.4441 










4 


-1-0 1 734 


—9 7498 


-H 1 909 


— 10 736 


-1-9 5963 


n 




5 


+0.1837 


-0.0875 


+0.0158 


-0.1350 










6 


-0.9877 


-0.8168 


+6.1555 


-4.4380 










7 


+0.9644 


+0.6888 


-5.6151 


+4.0284 










8 


+0.0529 


+0.6318 


-3.3995 


+2.5754 










9 


-0.0764 


+0.3017 


-2.9632 


+1.9311 










10 


+0.0690 


-1.5360 


+9.5808 


-10.168 


+2.9807 







11 


-1.4095 


-0.3543 


+6.4780 


-4.7956 










12 


-0.0043 


-0.3779 


+7.7358 


-42.390 


+61.401 


-26.778 
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FIGURES 

Light element abundances ^(77) ± 2(j ] (r)) 
1 1 1 1 1 — i — i — i — - 

' ~ 4 He moss fraction ~ 




FIG. 1. Primordial abundances (solid lines) and their ±2a bands (dashed lines), as functions 
of the nucleon-to-photon ratio rj. 
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FIG. 2. Logarithmic derivatives \it of the abundances Y{ with respect to the reaction rates Rk, 
as functions of 77. 
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Abundance uncertainties and correlations 

(2, 3, 4, 7 = D, 3 He, 4 He, 7 Li) 




FIG. 3. Fractional abundance uncertainties (JijYi (upper panel) and their correlations pij (lower 
panel), as functions of rj. 
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FIG. 4. Monte Carlo estimates of o-j/T, (SKM '93 |j, dots) and (KK '94 (2|, dots), 
compared with our analytic evaluation (solid lines), using the same inputs. 
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Standard BBN predictions vs data (1a errors) 




x = log 10 (7?X10 10 ) 

Y 2 = D/H 

Y 4 = 4 He mass fraction 
Y 7 = 7 Li/H 

Data set A 

Y 2 = (1.9±0.4)X1CT 4 
Y, = 0.234Y0.0054 
Y 7 = (1.6±0.36)X1CT 10 

Data set B 

Y 2 = (3.40±0.25)X1CT 
Y 4 = 0.243Y0.003 
Y 7 = (1.73±0.12)X10~ 




2.3 2.4 2.5 

10 x Y. 



1 x Y c 



FIG. 5. Standard BBN predictions (dotted lines) in the 2-dimensional planes denned by the 
abundances Y2, I4, and Y7, as functions of x = log 10 (??/10 _10 ). The theoretical uncertainties are 
depicted as la error ellipses at x = 0, 0.1, 0.2, . . . , 1. The crosses indicate the two observational 
data sets (with la errors). 
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Standard BBN fit to experimental data 




FIG. 6. Our x 2 fit to the data sets A and B, including observational and (correlated) theoretical 
errors. 
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Big Bang Nucleosynthesis — 
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FIG. 7. Individual contributions of different reaction rates Rk to the uncertainties in I2, ^4, 
and Y7, normalized to the corresponding total errors a 2, 04, and cr 7 , for rj = 1.78 x 1(T 10 . Each 
arrow corresponds to the shift SYi induced by a +lcr shift of Rk- Some small error components 
have not been plotted. 
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Big Bang Nucleosynthesis — Error Components 

at ri = 5.13X1CT 10 
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FIG. 8. Same as Fig. g, but for rj = 5.13 x 10" 
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